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Abstract 

We explore the dark matter sector in extensions of the Minimal Supersymmetric Standard Model 
(MSSM) that can provide a good fit to the PAMELA cosmic ray positron excess, while at the same 
time addressing the little hierarchy problem of the MSSM. Adding a singlet Higgs superfield, S, 
can account for the observed positron excess, as recently discussed in the literature, but we point 
out that it requires a fine-tuned choice for the parameters of the model. We find that including an 
additional singlet, ^, allows both a reduction of the weak-scale fine-tuning, and an interpretation 
of the cosmic ray observations in terms of dark matter annihilations in the galactic halo. Our setup 
contains a light axion, but does not require light CP-even scalars in the spectrum. 
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I. INTRODUCTION 



A. Dark Matter and Cosmology 

Within the standard A Cold Dark Matter cosmological model (ACDM), which repro- 
duces the available experimental observations with remarkable success, about a fifth of the 
energy density in the universe is contributed by dark matter (DM) pQ. The DM cannot 
be in the form of ordinary baryonic matter as deduced from considerations of cosmological 
nucleosynthesis and observations of the anisotropies of the cosmic microwave background 
(CMB). Despite constituting most of the mass in the universe and playing a crucial role 
in the growth and dynamics of structure, either the lack of appreciable interactions with 
observable particles or its mass have so far prevented the determination of its origin and 
composition. 

On the other hand, elementary particle theory provides several candidates for the DM, 
such as Weakly Interacting Massive Particles (WIMPs). What makes WIMPs attractive DM 
candidates is that their existence is motivated independently in particle physics extensions 
that address some of the shortcomings of the Standard Model. For instance, one of the 
by-products in extensions relying on low energy supersymmetry (Susy) to alleviate the 
hierarchy problem in the higgs sector of the SM, is the natural occurrence of a stable WIMP. 
The neutralino in the Minimal Superymmetric Standard Model (MSSM) is the simplest 
viable example. By virtue of its weak interactions, the thermal relic density of a WIMP is 
of the right order of magnitude, and its existence can be probed through both direct and 
indirect detection experiments (see e.g. (2J |3] for reviews). 

Although no convincing signal has yet been found, the latest PAMELA data [H E] , sug- 
gesting a new source of galactic positrons, has been interpreted as the result of annihilation 
of DM particles in the galactic halo [SHIS] or DM decay [MEo] . 

On the other hand, several astrophysical sources such as pulsars [26-f2"5] . supernova rem- 
nants [291 ED], or secondary production in regions where cosmic rays are accelerated [3TH31] 
may also account for part of, or perhaps all, of the observed fluxes. In addition, the hard 
injection spectrum required to fit the positron fraction measured by PAMELA, and the non- 
observation of an equivalent anti-proton excess [35] exclude annihilation of thermal WIMPs, 
such as the neutralino in the MSSM, as a viable explanation. 
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A DM interpretation of the data demands that the particles making up the dark galactic 
halo annihilate mostly to charged leptons with a cross-section which is ~ 10 — 100 times 
larger than the canonical value leading to the correct cosmological abundance via thermal 
decoupling [36J, (crv) ann ^3x 10 -26 cm 3 s. MSSM neutralinos could annihilate primarily to 
leptons, either due to radiative corrections in a small region of parameter space [37J, or for a 
wino-like LSP of about 200 GeV as discussed in [38], although a non-standard cosmological 
evolution has to be invoked to explain the boosted annihilation cross-section required by 
observations. A nearby clump of DM could raise the annihilation rate [3D], but it is highly 
unlikely that a sufficiently large clump can be found in the solar neighborhood [ID]. The 
generic expectation in the MSSM is that neutralinos annihilate to a mixture of heavy quarks 
and Higgs bosons producing a spectrum that is too soft to account for the PAMELA data, 
while also over-producing anti-protons |41j . 

Several alternative particle physics scenarios have been proposed that predict DM parti- 
cles can explain the rise in the positron flux at high energies without conflicting with other 
measurements. For instance, if the DM particles annihilate with a weak-scale strength to 
light metastable mediators which are very weakly coupled to the Standard Model, then 
kinematical constraints preclude any final states other than light leptons. In addition, the 
exchange of light mediators results in long-range interactions that enhance the annihilation 
cross-section when the DM particles are moving at non-relativistic velocities, as is the case 
for present day processes occurring in the galactic halo. This so-called Sommerfeld enhance- 
ment [DJ, H21 E3], reconciles the required underlying weak-scale interaction at the time of 
DM freeze-out that generates the correct relic abundance with a much larger annihilation 
cross-section in the present galactic environment. The scenarios in [Dt lT2"lH3] . among others, 
provide particular implementations of these principles, but they are not primarily motivated 
by a solution of the hierarchy problem. 

B. Natural Models of Dark Matter 

In this paper, we seek a model of dark matter that addresses the hierarchy problem and 
other naturalness constraints while still generating the fluxes observed by PAMELA. We 
require the following 

• Stable dark matter with a mass of 100 GeV or above. 
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• Dark matter annihilations dominantly to light leptons. 

• Solution to the hierarchy problem. 

• Minimal fine-tuning among model parameters. 

• Thermally generated dark matter. 

• Passes constraints from accelerators. 

We find that this set of requirements will place severe restrictions on the form of the model. 

We begin with a consideration of naturalness in supersymmetric models. The MSSM, 
while solving some problems of the Standard Model, introduces new theoretical problems 
that suggest it might not be a complete description of physics at the electroweak scale. It 
must contain a mass term \x for the two Higgs doublets, H u and H^, which can neither vanish 
nor be naturally large (~ Mqut or ~ Mpi) for phenomenological reasons. The lack of an 
explanation of /j, ~ Ms usy constitutes the /i-problem of the MSSM [45J. The addition of a 
singlet chiral superfield, S, to the particle content of the MSSM provides an elegant solution 
to the /i-problem: the scalar component of S obtains a vacuum expectation value (vev) of 
the right order, dynamically generating the mass term fi. 

Another challenge to these models is that the non-observation of the Higgs boson at 
LEP-II requires large soft-supersymmetry-breaking mass parameters to raise the mass of 
the lightest CP-even Higgs above the tree-level prediction, rrih < mz, in the MSSM. The 
discrepancy between the large size of these soft Susy breaking terms compared to their 
natural scale, the electroweak scale, is known as the little-hierarchy problem. Here again, 
the scalar components of S, that mix with the neutral scalar components of H u and Hj, 
can alleviate the little fine-tuning problem of the MSSM |47)H4*5] by lifting the Higgs mass or 
allowing for new Higgs decay modes that weaken the LEP-II limits [48-50J. The resulting 
model is the Next-to- Minimal Supersymmetric Standard Model (NMSSM, see e.g. [511 E2] 
for recent reviews), and its dark sector can differ considerably from that of the MSSM. Much 
like in the Higgs sector, mixings with the fermionic component of S, the singlino, result in an 
extended neutralino sector. The lightest neutralino can have a sizeable singlino component 
and, if it is the LSP, the expected signatures at colliders and the DM phenomenology could 
markedly differ from the minimal scenario. 
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Interestingly, it has been pointed out that the positron excess observed by PAMELA 
can be explained by neutralino annihilation in the NMSSM [ToT - [T%] . The richer Higgs and 
neutralino sectors can potentially accommodate the ingredients shown in [HI [T21 HI] to result 
in enhanced mostly leptonic fluxes. The fact that this scenario is motivated independently 
from particle physics considerations as outlined above, makes it even more appealing. It is 
this last point that we set forth to study in this paper. We revisit in Section [IT] the region 
in the NMSSM parameter space that allows for an explanation of the reported cosmic-ray 
anomalies in terms of neutralino annihilations, and study whether the original motivation 
of naturally addressing the little higgs and /^-problem is preserved. Our findings show that 
this is not generically possible without accidental relations among the parameters, losing 



the naturalness motivation. In Section III we add another singlet superfield, to the dark 
sector, which suffices to avoid reintroducing fine-tunings in the electroweak sector. A study 
of a similar model, with vector-like dark matter, was presented in [12] in the context of 
an axionic sector. Here we do not require a light singlet scalar to generate a Sommerfeld 
enhancement, as we take a light pseudoscalar to be sufficient. We study the dark matter 
sector in Section |V] and we find somewhat different behavior in our scenario than the one 
presented in [12] regarding the behavior of the extended NMSSM model. Our results are 
summarized in Section IVTl 



II. DARK MATTER PHENOMENOLOGY IN SINGLET EXTENSIONS OF THE 
MSSM 

In supersymmetric extensions of the SM the masses of the up-type quarks and down-type 
quarks are generated by the vevs of two Higgs SU(2)-doublets H u and H^- With this minimal 
field content in the Higgs sector, as found in the MSSM, a dimension-full coupling fiH u H v 
must be appear in the superpotential. For phenomenological reasons the parameter /i has to 
be of the order of the electroweak breaking scale, t>Ew, which is orders of magnitude below 
the natural value in the MSSM, [i ~ A, where A represents the ultra-violet (UV) cut-off of 
the theory (GUT or Planck scale). In the NMSSM (see e.g. [SH E2] for recent reviews and 
further references), the mass term /i is replaced by the vev of a scalar field - induced by 
the soft Susy breaking terms - which has a Yukawa coupling to the Higgs doublets. The 
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simplest scenario allows only scale-invariant and renormalizable terms in the superpotential: 

W = \SH u H d +^S 3 , (1) 

where A, k are dimensionless, and we have not included couplings to the lepton and quark 
matter fields of the MSSM, which are not relevant here. Electroweak symmetry breaking 
generates vevs for the scalar components of both MSSM Higgs superfields, H u and H4, as 
well as a vev (S) for S of the order of the weak scale. In this way the NMSSM generates an 
effective \i paramater, 

li = X(S), (2) 

of the order of the weak scale, t>Ew ~ 174 GeV, which solves the /i-problem of the MSSM. 
The scalar components of S mix with the neutral scalar components of H u and H^, and this 
results (in the absence of explicit CP violation) in three CP-even and two CP-odd neutral 
scalars. Mixing of the fermionic components of S with the neutral higgsinos and gauginos 
leads to five neutralinos in the spectrum of the NMSSM. The tree level Higgs and neutralino 



mass matrices can be found in e.g. J5TJ [52], or deduced from the expressions in Section III 
by taking the limit £, — >■ 0. 

As mentioned above, the richer Higgs sector in the NMSSM allows for a resolution of 
the little fine-tuning problem of the MSSM [45H5TJ1 1515] . In addition, the appearance of new 
Higgs decay modes, together with a possible singlino component, can result in a markedly 
different phenomenology for the light neutralino, which is a candidate for DM as in the 
MSSM. Interestingly, it has been pointed out by several groups that an explanation of 
the anomalous cosmic-ray fluxes reported by PAMELA in terms of NMSSM neutralino 
annihilations is possible. 

In the scenario envisaged in [16], the LSP is a neutralino of the bino-type with Higgsino 
mixings. This allows for a large annihilation cross-section into the lightest CP-even plus CP- 
odd scalars, Xi + Xi ~ > ^1 + As explained below, it is technically natural for a CP-odd 
scalar in the NMSSM to be light, and this restricts the subsequent decays to a\ — > ji + pT . 
Hence, a flux of energetic positrons is obtained, while kinematics prevents the production 
of anti-protons, in agreement with observations. The positron fluxes can only be brought 
in agreement with the data, if the mass of the heavier CP-odd scalar, 02, is dialed to be 
m a2 w 2m x o. Then, neutralino annihilation proceeds through the resonant 02 state, which 



would also increase the intensity of monochromatic 7-ray lines [51] • This scenario exemplifies 
the possibilities allowed by the larger Higgs sector of the NMSSM, but the careful choice of 
parameters required for the resonant annihilations to occur opposes the initial motivation 
to go beyond the MSSM. 

These models introduce further theoretical questions. For instance, the annihilation 
cross-section is more than two orders of magnitude larger at all times, so that a non-thermal 
production mechanism has to be invoked to match the observed DM density Since a non- 
thermal neutralino in the MSSM is also capable of explaining the PAMELA signal [3E], there 
would be little motivation to go beyond the MSSM. Furthermore, in the scenario presented 
in [16] , the loop corrections to the mass of a\ are not suppressed by any small parameter, so 
a fine-tuning is required to achieve the sub-GeV mass in their parameter set that primarily 
produces muons in dark matter annihilation. The neutralino is also required to be less than 
the top mass, which may present a difficulty in fitting the shape of the PAMELA spectrum. 

A different possibility, considered in [17], is that the neutralino could have a sizeable 
singlino component in the NMSSM. Annihilations would still proceed to the lightest scalars 
h\ + oi, both of which would now be mostly singlet-like, not just a%. Due to its very small 
couplings to all quarks, leptons and gauge bosons, a singlet-like h\ is not bound by the LEP- 
II limit, rnu SM > 114 GeV, and could in fact be much lighter. Such a light scalar would 
naturally enhance the rate of neutralino annihilations in the Galaxy through the Sommerfeld 
enhancement, without the need of resonant annihilations and a non-standard cosmo logical 
evolution. Hence, in this setting all the ingredients of the DM scenarios engineered to 
reproduce the PAMELA signal j9j H21 01] would be realized in a framework motivated 
independently from considerations of naturalness of the electroweak interactions. One should 
still make sure, however, that this does not require particular choices of the parameters to 
solve the /i-problem, which was the main motivation to go beyond the MSSM [65J. 

For a sufficient enhancement of the DM annihilation cross-section, the mass of the lightest 
CP-even scalar, h\, can be estimated as m/ ll < At 2 m x o/47r, where At is the trilinear coupling 
in Eq. ([I]). The anomalous trend in the positron fraction reported by PAMELA extends 
up to energies of ~ 100 GeV, which also sets the (minimum) mass of the neutralino. For 
moderate values of At < 1 the mass of hi should be roughly m hl < 10 GeV (more details can 
be found in P EEZ])- Since the decay of the lightest CP-odd scalar, a 1; should yield mostly 
charged leptons (e^ 1 or /r* 1 ), we demand m ai < 1 GeV. We also require that the /1 term is 
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dynamically generated as in Eq. ([2]). 

The parameter space of the NMSSM defined by Eq. (IT]) can be studied with the package 
NMSSMTools [SS], which calls MicrOMEGAS [67] to calculate the relic density [22 EE]. 
Given a point in the NMSSM parameter space, the program suites check its viability against 
negative particle searches and rare-decay bounds from accelerators. We performed a number 
of scans, both over random points and on grids, in different regions of parameter space. For 
instance, in the region defined above, which is the region considered by [16], we randomly 
selected one million points that had A < 0.6, with a flat distribution, and k < A. We fixed 
the soft terms A\ and A K to have an exponential distribution peaked at zero, with a width 
of 10 GeV, to attempt to take advantage of a U(1)r symmetry that makes a pseudoscalar 
light, as in [16]. We also fixed tan/3 between 2 and 10, and fi between 100 and 240 GeV, and 
we randomly assigned all free signs to be positive or negative. The results demonstrate the 
difficulty in finding valid models. Approximately 0.2% of the parameters passed all present 
experimental constraints. Of these, none had a pseudoscalar below 1 GeV. Additionally, the 
neutralino had a mass of order the light scalars, 0(10 GeV) rather than the required 100 
GeV or greater. The spectrum of points that passed LEP and other constraints is shown in 
Figure [TJ 

Random scan of points (pass constraints) 




10 20 30 40 

Mass of a ! (GeV) 

FIG. 1: Values of the mass of the light pseudoscalar and the neutralino for parameter sets that 
pass constraints for the scan described in the text. 

On the other hand, non-scale invariant terms in the superpotential were allowed in [T7] . 
Including supersymmetric mass and tadpole terms, fiH u Hd + ^fi'S 2 + £fS, it is possible 



to viable points in this extended parameter space [18]. Even though, some of these terms 
are non-vanishing in various scenarios, they contradict the theoretical motivation for the 
NMSSM that the /i-term should be dynamically generated as in Eq. (|2]). We cannot exclude 
the possibility that a more throughout scan might reveal a valid region that meets our 
stricter requirements, but the general considerations below should illustrate the difficulties 
for a generic natural extension of the MSSM. 

Light CP-odd scalars can easily appear in the NMSSM in the form of (pseudo-)Nambu- 
Goldstone bosons associated to approximate global symmetries. Both, in the R-symmetry 
limit, A\,A K — > 0, considered in [TEHTB"] . and in the Peccei-Quinn (PQ) limit, k — > 0, 
on which we focus below, axion-like particles would be naturally present in the spectrum. 
Concrete scenarios exist (as reviewed in e.g. [5T|l52]). in which this symmetries are explicitly 
broken, yet the axion remains light. The anti-proton constraints can, thus, be naturally met 
in the NMSSM. 

There is no symmetry, however, that prevents the mass of the CP-even scalar, hi, from 
growing above the scale required for the Sommerfeld enhancement to be effective. In fact, 
as shown in Section III, requiring that the /i-parameter is generated as in Eq. g) and that 
the lightest neutralino be heavier than ~ 100 GeV, cannot be simultaneously achieved while 
keeping m/ ll < 10 GeV. If we allow for non-scale invariant terms in the superpotential as 
in [TTJ [TS], the lightest scalar is no longer linked to the electroweak-scale /i-term and can 
be lighter. However, this is not a natural outcome of singlet extensions of the MSSM, and 
requires the tuning of the parameters. 

The main reason to include a light CP-even scalar is to enhance the annihilation cross- 
section in the low velocity limit, so that the positron fluxes from annihilations in the Galaxy 
can be large enough, while obtaining the right relic density from thermal freeze-out at 
the time of decoupling. As outlined in [91 H2J H3], such a light particle mediates long 
range forces that generate the non-perturbative Sommerfeld enhancement. A light vector 
particle also results in long-range forces, and this was in fact the first scenario in which 
the importance of this effect for DM was studied [12]. Other possibilities, such as CP-odd 
mediator particles, have been dismissed, since they do not give rise to spin-independent long 
range interactions at tree level. However, since pseudo-scalars do give rise to long-range 
forces at the one loop level [57], this conclusion should not necessarily apply in the non- 
perturbative regime. Indeed, as shown in [58] . (pseudo-)Nambu-Goldstone bosons can cause 
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a Sommerfeld enhancement. As discussed above, we do have light pseudo-scalars in our 
spectrum, and we might ask whether lifting the limit on the mass of the hi, allows a DM 
explanation to the PAMELA observations, while abiding to Eq. (|2]). 

Unfortunately, we find that this is not sufficient. In order to generate a splitting between 
the scalar and fermion components of the S superfield, in general we need to increase the 
mixing with the component of H u and Hd by increasing A. However, this also increase the 
magnitude of radiative corrections to the mass of the light scalars, requiring a high-degree 
of fine-tuning to maintain a sub-GeV pseudoscalar mass. One might hope to escape this 
problem by introducing additional superpotential terms, as in [17] . however a subsequent 
study of the parameter space by (H] found that this scenario requires A ~ 1CT 3 , a small 
coupling with no clear origin. 

The previous discussion shows the difficulties in constructing a DM sector that allows an 
explanation of the rising positron fraction, and that is motivated by naturalness considera- 
tions of the electroweak interactions. A possible solution, with only scale invariant terms, 
was put forward in [12] , where it was noted that an additional singlet could be added to the 
PQ-limit of the NMSSM. In the following we consider an scenario in the spirit of p2], and 
study the restrictions imposed by Eq. (J2|. In addition, we note that a similar superpotential 
was proposed in [59] to extend the NMSSM by a right-handed neutrino superfield, which 
serves as a suitable thermal dark matter candidate. As the motivation in [59] is not to ex- 
plain observations made by PAMELA, the authors examine a different region of parameter 
space and an additional superfield coupling that is not relevant for our purposes. 

III. ADDING A DARK MATTER SINGLET 

Following from the above consideration of the NMSSM, we argue that a super symmet- 
ric model that is not fine-tuned and contains leptophilic dark matter annihilation requires 
extending the MSSM by more than one field. Here we will consider the simplest case: the 
MSSM with an additional two superfields: a gauge singlet S, as in the NMSSM, and an 
additional singlet ^ that will serve as dark matter. In our arrangement, the S and ^ to- 
gether form a dark sector that is moderately secluded from the MSSM fields. While at first 
glance this theory may appear to be overly general, we will show below that observations 
and naturalness significantly restrict the allowed couplings. 
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To naturally generate a light pseudoscalar, we will take the U(1)pq limit of the NMSSM, 
but will include an explicit symmetry breaking term k. We take this limit both because it 
generates a light pseudoscalar with a controllable mass, and because that small mass is more 
stable against radiative corrections than the one provided by U(1)r symmetry. We assume 
the \& superfield has a Z2 symmetry which will guarantee a stability 

The superpotential of this theory is 

W =^ 2 + XSH u H d +^S 3 , (3) 

where S is the singlet of the NMSSM, and H u and Hd are electroweak Higgs doublets, as 
in the MSSM. The new superfield 4/ consists of a fermionic component ip, which will be a 
stable dark matter candidate, and complex scalar <fi. We have not included the other matter 
fields of the MSSM, which are not relevant here. The corresponding scalar soft terms are 

Koft = ^S0 2 - \A x SH° u H° d + ^A K S 3 + h.c, (4) 

where we have dropped electrically charged components. Since k is a small breaking of 
the U(1)pq symmetry, we will take /«< 1. The typical sizes of the other parameters are 
discussed below. 

As in the NMSSM, we require that the effective /1 parameter given in ^ is of order 
the electroweak scale. Electroweak symmetry breaking does not generate a VEV for the 
scalar partner of the dark matter, <fi. The consistency of this minimum can be checked by 
examining the eigenvalues of the Hessian for the scalar potential. In this paper we will not 
study the global structure of minima. 

IV. MASS SPECTRUM 

Having defined the theory, we will turn our attention to the mass spectrum of the theory 
in the presence of the additional states. We will denote the real scalar and pseudoscalar 
components of the </> and S fields with the subscripts s and a, respectively, and use the usual 
notation for the Higgs VEVs, 

(H° u ) = v u , (H° d ) = v d , (5) 
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where v 2 + v d — i>| w . In reference to the MSSM, it will be convenient to define an effective 
B parameter, 



M 2 RS 



B = A x + k(S) = A x + k^. 

A 

The CP-even mass matrix, in the basis (H d , H u , S s , (fi s ), is given by 

( g 2 v 2 d + Bfi^ v u v d (2X 2 - g 2 ) - Bfi 2\{iv d - v u (XB + K fi) 
g 2 v 2 u + BpL^ 2\/iv u - v d (XB + K fi) 



(6) 



v„v d (±f - k\) + kA k * 









V 



(7) 



where we drop terms of 0(k 2 ). 

The Higgs sector contains a neutral massless Goldstone mode, G, that becomes the 
longitudinal mode of the Z boson. This can be rotated away, leaving a single physical Higgs 
pseudoscalar from the MSSM, A. After performing this transformation, the CP-odd mass 
matrix, in the basis (A, S a , a ), is given by 



Mf s = 



( MMS2^±^l (5 A - 3^ cos + sin /3) 

^(5A + 3^)-3fA^ 








\ 



V 



(8) 



j2^uV d \ 3 - /i(AA ? + Kfl- J 

The neutralino mass matrix is extended from the MSSM by the fermionic components of 



the superfields S and ^. In the basis (B, W°, H d , H u , S, ip), it is given by 

f M ± -^f 9 -^f \ 



M- 



92V d 





-/i -Xv u 
-\v d 
2f/i 



V 



(9) 



We will label the resulting scalar and pseudoscalar mass eigenstates with hi and Oj re- 
spectively, ordered from lightest to heaviest. Similarly, we will label the neutralinos with 
iVj. To avoid confusion, we will continue to label the states of \& with ip and 0, since these 
states are the same as their flavor states. 
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We will be interested in the region of parameter space where A is small, but not unnatu- 
rally so. In the PQ limit, which we consider, k is expected to be very small, but A is not. In 
this case, the lightest scalar and pseudoscalar, hi and ai, are composed primarily of singlet 
S states. The Higgs with significant couplings to the Standard Model are h,2, h% and ei2. 
The lightest neutralino, N\, is composed primarily of singlino S. Since the lightest states 
are all only slightly mixed with other Higgs, their coupling are very nearly those of the S 
flavor states. 

The behavior of dark matter in this theory primarily depends on the value of the masses 
of the light particles. In the region of parameter space we will look at below, the Sommerfeld 
enhancement is generated primarily by a\. The lightness of this particle also kinematically 
enforces the creation of light leptons rather than baryons in dark matter annihilation. To 
leading order in k and A, the tree-level mass of this particle is 

9 

m 2 ai = -vl w n\sm2(3. (10) 

The lightest scalar, h±, has a leading mass that is independent of k and so is generi- 
cally unsuppressed relative to m ai . The more complicated mixing of the 3x3 CP-even mass 
submatrix does not lend itself to compact expression for this mass. To leading order in A, 
the mass is given by 



rn 



(8A A /i 3 - 8A 2 /i 2 sin 2/3 + 2A X ^(A 2 X - 2g 2 v 2 ) sin 2 2/3 



2 -A 2 sec 2 2/3 
^ ~ 2g*A x n 

+ g 2 v 2 (Al + 4/i 2 ) sin 3 2/3) . (11) 

The hierarchy of scalar and pseudoscalar masses is illustrated for typical values in Figure [2j 
Positivity requirements for the mass place restrictions on the allowed values of the other 
coupling constant and soft terms. This relationship further implies the existence of a maxi- 
mal value for once the dimensionless coupling constants are fixed. An example of this 
for typical values is given in Figure [3j A \ is of the order of the typical soft breaking scale. 
A K is less relevant, as its contribution is suppressed by k. 

The scalar and pseudoscalar masses are shifted by loop corrections from the tree-level 
values given above. These corrections are suppressed by additional powers of A and k, plus 
loop factors. For A < 0.1, these corrections are sufficiently small that they do not affect 
the results we present. In particular, the light pseudoscalar receives mass corrections that 
are less than a GeV. Note that this is true in the PQ-limit with small k, but requires an 
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FIG. 2: Contour plots of the tree- level values of m 2 ^ and m| as a function of k and A for the case 
of n = 130 GeV 



m hl 2 (GeV 2 ) 




A x (GeV) 



FIG. 3: Values for the tree-level square mass of the light scalar, m? , where // = i>ew 5 A = 0.1, 
tan/3 = 4 and A\ = 500 GeV. 

unnaturally small value of A (which is not related to any enhanced symmetry) in the regions 
considered in [T6HT8]. 

The supersymmetric partner to these light scalars, N\, typically has a mass of a few GeV 
that is intermediate between hi and a±. To leading order in A, gi and g2, an approximation 
good to 10% or better, this particle has a tree-level mass 

AV (4M 1 M 2 /isin2/3 - 2(g^M 1 + g\M 2 )v 2 cos 2 2/3) 



mi 



(12) 
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Finally, there is the mass of the dark matter itself. The fermion ip is stable and has a 
mass that is fixed by (S) and the coupling £, 

m^ = £j. (13) 

Note that if \i is near the electroweak scale, then is 0(TeV) for the relevant values of 
A. We then get a dark matter mass of (9(100 GeV) or above for £ > 0.1. 



V. DARK MATTER 



Within the model described above, we identify several possible dark matter scenarios. 
First, we assume that supersymmetry breaking is gauge-mediated, so the true LSP of the 
theory is the gravitino. The stability of ip is protected by an unbroken Z2 symmetry, which 
makes it a suitable heavy dark matter candidate. The S and a have mass that is typically 
somewhat larger, so can decay through the channel Si(l — > ipNi. The N±, the NLSP of the 
theory, will then decay to a gravitino and a singlet scalar. 

Even if < + , the can still decay to a gravitino and a tp, either in tree-level 
processes with a virtual Ni, or in 1-loop processes with an h\ or an a\ in the loop. To avoid 
problems in the early universe, we assume either that is sufficiently massive to decay to ip 
and Ni on-shell, or that the suppressed processes allow the to decay prior to the BBN. 

If the are long-lived due to suppressed decay processes, a careful consideration of the 
early universe is required. If a significant number of particles remain after freeze-out, a 
late decay to gravitinos after nucleosynthesis will ruin the successful predictions of early 
universe cosmology. However, the has additional self-annihilation channels generated by 
3- and 4-point scalar interactions that tend to raise its cross-section significantly, and their 
relic density becomes irrelevant. If the mass is nearly degenerate with the ip mass, then 
there may be additional cross-annihilation channels that will shift the ip relic density from 
the values given below. This offers the interesting possibility of breaking the relationship of 
m^p and the relic density, both of which are set by £. 

There is also the possibility that supersymmetry breaking is gravity-mediated. In this 
case, the singlino is the true LSP and absolutely stable. This is problematic in the model 
presented above, as the dominant channel for singlino self annihilation, to Higgs through an 
s-channel, scales as k 2 , and so JVi is overproduced thermally in the early universe. However, 
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in a slightly expanded version of the model, in which we allow dimensional superpotential 
terms, we can arrange for a U(1)pq symmetry with non-infinitesimal k. A version of this 
model was described in [T7] . Here the self-coupling can be much larger, and we can arrange 
for a two-component dark matter system that is a mixture of the heavy if) and the light N\. 
Since this scenario requires an extended parameter set, we will not consider it further here. 

We will make the new fields moderately secluded through a small, but natural, value of 
the coupling A ~ 0.1. Above this value, the radiative corrections to the scalar masses shifts 
their values from the tree-level significantly. Below this value, the annihilation cross-section 
is suppressed, and the universe tends to overdose. After A and fi are fixed, the coupling 
£ is still free to fix both the dark matter mass and its annihilation cross-section. We will 
show below that £ ~ 0.4 is an appropriate value. Our couplings are arranged in the loose 
hierarchy 

k < A < f . (14) 



A. Relic Density 

The relic dark matter is thermally produced in this scenario. While the coupling of 
the dark matter to the standard model fields is small, the coupling £ of the dark matter 
particles to fields of S is essentially unconstrained. The value of this parameter fixes the 
annihilation cross-section, and it is notable that the same values that produce dark matter 



mass of (9(100 GeV) in (13) also yield the observed relic density. 

A number of annihilation channels are potentially relevant. There is an s-channel process 
mediated by hi or a\ with a two-particle final state consisting of a combination of scalars 
hi, di, h 2 or h s . Higgsino production in this channel is possible, but suppressed by k if it is 
allowed kinematically. 

There are two relevant t- and u-channel processes. The first is mediated by and results 
in two Ni particles in the final state. The second is mediated by if) and results in an hi and 
ai. The processes that are not A-suppressed are shown in Figure |4j 

When comparing the cross-sections of the various channels, one might expect that the 
Higgs production is suppressed by a large amount, as these matrix elements appear to be 
proportional to A 2 . However, some channels are enhanced by (S), which is proportional to 
A -1 , and they are only suppressed by one power of A. Since we do not have £ ^> A, we can 
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{h 1 ,a 1 ,a 1 } 



{hi,ai,hi} 




/ii,ai \ 



{hi,hi,ai,ai,h 2 ,h 2 ,h 3 } 



{h 2 ,h 3 ,h 2 ,h 3 ,h 3 ,h 2 ,h 3 } 



FIG. 4: Processes that contribute to the annihilation of tp which are unsuppressed by k. 



not a priori ignore Higgs production during ip annihilation. 

We implemented this model and evaluated cross-sections within CalcHEP jHO] and found 
the annihilations are dominated by the process ipip — > NxN\. This non-relativistic cross- 
section, as a function of \x and ip, is shown in Figure [5} The process ipip — > h%a\ provides a 
significant subdominant contribution, approximately 10% of the previous process. A plot of 
the cross-section for the same region of parameter space is given in Figure |6j 

Higgs production s-channel processes to aih 2 and dih 3 are typically suppressed by another 
factor of 10 2 from h\a\ production, and so can be safely ignored for the parameter space of 
interest. 

Finally, we note that in this model the secluded sector may decouple from the other 
fields of the NMSSM and undergo a separate evolution in the early universe. If there are 
additional heavy degrees of freedom in the secluded sector that come out of equilibrium 
during this period, the temperature of the secluded sector will diverge from the temperature 
of the photons. This shift in temperature may alter the relationship of cross-section and relic 
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FIG. 5: av for the annihilation of tptp to N\N\ in units of cm 3 /s with A = 0.1. The gray region 
shows the cross-section that generates the observed relic abundance. 
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FIG. 6: <tk for the annihilation of ?/;^ to /iiai in units of cm 3 /s with A = 0.1. 

density through an effect like the one described in [TT]. Here ip and a\ would, respectively, 
play the role of X and Y from that paper. The Zi is provided by <p a and <f> a , and the 
temperature difference could be enhanced by including additional heavy unstable singlets. 
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B. Galactic Annihilations 



In the galaxy, the effective cross-section for dark matter annihilation is increased due to 
the Sommerfeld effect. In our model, this effect arises predominantly from the exchange of 
the pseudoscalar a%, as described in [58] . 

As described in the preceding section, dark matter annihilation will produce S particles 
that quickly decay to oi, or will produce a\ and hi directly. The zero-temperature cross- 
sections, relevant in the galaxy, for dark matter annihilation to these particles have fairly 
simple approximate analytic expressions, 



£ 2 A 2 

<™W->M = I^^ + 0(A4) ' (lf)) 

i2\2 

2048vr/i 2 

When making numerical calculations, we do not use these approximate expressions, but 
instead keep all orders of A and T. 

The h\ and a\ produced in these annihilations have small mixings with the MSSM Higgs 



av^^ hiai = —— — - + 0(A 4 ). (16) 



that will cause them to decay to the heaviest allowed states. The mass of cti, given (10), can 
be kept sub-GeV for sufficiently small k, and loop corrections to the mass can be kept under 
control with a moderately small A. This allows kinematic enforcement for the production of 
light leptons during annihilation, as in the mechanism of [9j [11] and others. This can give 
rise to a flux of positrons that is measurable about the modeled astrophysical background, 
as was observed in the PAMELA experiment [3]. 

The hi typically has a mass above the threshold for baryon production, and so will con- 
tribute to the anti-proton flux. However, the cross-section of this subdominant annihilation 
to baryons is sufficiently small that the excess of anti-protons will lie well within the current 
experimental uncertainty. 

We estimate the positron spectrum within the energy range measured by PAMELA by 
using approximate solutions to the diffusion equation, 

/ - K(E) " V 2 / - ^ (b(E)f) = Q, (17) 

where / is the positron distribution, K is a diffusion constant, b is an energy loss coefficient, 
and Q is the positron source term. For the case of cylindrical symmetry, an analytical 
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solution for the positron flux observed at Earth is given by 

- B sSm ^)lz dEav 1it I(XD{E ' E)) - (18) 

where K is the diffusion coefficient, p & is the dark matter density at our solar system, / is 
the halo profile, and A^> is the positron diffusion length. An overview of the approximate 
solution to this equation in the presence of cylindrical symmetry is given in [61 J, and we use 
the "medium" NFW dark matter profile from that paper. 

The electrons and positrons are created at the end of a chain decay. In the dominant 
process, the if) annihilation produces two N\ particles, each of which decays to a gravitino 
and an a%. Depending on its mass, the a\ decays to muons or electrons. We numerically 
calculate the resulting spectrum, dN e + / dE e + by selecting random directions for each step 
of the decay, isotropic in the decaying particle's rest frame, and then boost the electron 
4-momentum to the galactic frame. We perform 10 5 trails, and fit the resulting spectrum to 
to the three-parameter function C\E C \ + C3, which yields excellent agreement. We perform 
this procedure for each combination of fi, A and £ that are used below. 

For each point in parameter space, we fit the positron spectrum generated by dark matter 
annihilations to the PAMELA spectrum with an overall scaling "boost factor." The observed 
spectrum, originally reported in jl], has recently been updated using a new statistical method 
[5]. The method has resulted in a spectrum that is slightly softer, which suggests a larger 
dark matter mass and reduced boost factors. We fit our model to both versions of the data. 
A typical spectrum that results from the numerical fitting procedure is shown in Figure [7} 
The boost factors required may be a combination of enhancement in the local annihilation 
due to structure within the dark matter halo, and the pseudoscalar Sommerfeld effect. A 
plot of the required boost factors across a portion of parameter space is given in Figure [8} 
The large hierarchy between pseudoscalar mass and dark matter mass, rriff /m ai ~ 0(500) 
easily provides the boost factors required for jj, < 100 GeV. 

VI. DISCUSSION AND CONCLUSIONS 

Our exploration of the dark matter sector in singlet, S, extensions of the MSSM shows 
that explaining the anomalous positron fluxes observed by PAMELA with neutralino anni- 
hilations is not possible without a carefully tuned choice of parameters, in regions where the 



20 



Positron Ratio 
0.20 r 



0.15 



0.10 



0.05 




20 



40 



60 



80 



100 



120 



E (GeV) 



FIG. 7: The positron spectrum observed near Earth, from the reanalyzed PAMELA data given 
in [5], shown as a ratio of positrons to the sum of positrons and electrons. The points are the 
spectrum observed by PAMELA, and the line is the spectrum generated by the present model for 
the particular choice /u = 70 GeV, £ = 0.3 and A = 0.1 (which yields = 210 GeV), and a boost 
factor of 42. 
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FIG. 8: The boost factor required, as a function of £, to fit the model to the PAMELA spectrum 
for A = 0.1. The solid curves result from a numerical fit to the newer PAMELA [5], and the 
dashed curves result from a numerical fit to the original PAMELA spectrum [5j . Each neighboring 
pair of solid and dashed curves is for a different value of /i. From bottom to top, these are 
H = 70,90,110 GeV. 
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original particle physics motivation for these extensions is lost. 

In the NMSSM, the simplest such extension that addresses the little-fine tuning and the 
\i problems, it is hard to reconcile a neutralino at the electroweak scale, > 100 GeV, 
with the light (pseudo-)scalars that are required to avoid large anti-proton signals, and to 
enhance the present day annihilation rate. Indeed, this can only be achieved by adding non- 
scale invariant terms to the superpotential. These terms could emerge from an unspecified 
higher energy theory, but then the scale of the \x term cannot be explained in the framework 
of the NMSSM. 

We have showed how adding a dark matter singlet superfield, ^, to the NMSSM in the 
PQ limit, makes it possible to obtain a large positron flux from DM annihilations, without 
an unnatural choice of parameters. The DM, ip, annihilates primarily into a pair of NMSSM 
singlino-type neutralinos. These in turn produce pairs CP-odd higgs, which can be very 
light in the PQ-limit and would mostly generate light leptons upon decay. Anti-protons in 
this scenario are only generated by an annihilation to light scalars that include the hi, not 
in the annihilation to neutralinos. Hence the lack of antiprotons in the spectrum is naturally 
generated. 

A similar scenario had been considered in [12], but our setup differs in a few important 
aspects. Our DM is a Majorana fermion, whose stability is ensured by a Z2 symmetry, that 
annihilates primarily into a pair of NMSSM singlinos, whereas vector-like DM annihilation 
can be dominantly to hi + cti [12] • Positrons are generated by the decay of the a 1( which ap- 
pears further down in the reaction chain in our setup, and resulting in a softer spectrum. The 
updated positron spectrum is softer at high energies [5] than suggested by initial estimates, 
and we have shown how our model can fit the data with natural values of the parameters in 
the DM sector, £ < 1, and in the electroweak sector, A < 1 and fi ~ 100 GeV. The excess 
in the positron fraction could extend to higher energies, and there have been conflicting 
measurements of a large bump by ATIC [62J, which is not seen by Fermi [03] • Systematic 
uncertainties make these measurements challenging, and more experimental work needs to 
be done to clarify the spectrum above the PAMELA cut-off [51] . In view of these facts, 
we have not attempted to fit any of the higher energy experiments in our model. Let us 
briefly mention, that we can allow masses for the DM particle that are large enough to create 
positron fluxes at much larger energies. This is not the case in other scenarios, where going 
beyond the top quark mass is problematic [IE]. We also predict a softer spectrum than in 
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models where the positrons appear earlier in the annihilation chain [9j [121 HI fTTj . 

The annihilation rate in the galactic halo should be larger than at the time of decou- 
pling, and previous proposals have frequently invoked a light scalar particle to achieve this 
boost [91 H2J [T7J HI]. However, we have explicitly shown in a particular well- motivated 
setup that CP-even scalars are not expected to be light enough to generate the required 
Sommerfeld enhancement. Unlike previous studies, we do not have any in the spectrum, 
but we have pointed out that the light CP-odd higgs can also be used for this purpose. 
Hence, the relic DM is thermally produced in this scenario. While the theory may appear 
overly general, observations and naturalness considerations restrict the allowed couplings to 
a narrow range. The coupling £ of the DM particles to the singlet field S determines both 
the annihilation cross-section and the mass of the DM particles. It is remarkable that a 
value £ ~ 0.3 produces the required (9(100 GeV) mass, while it also yields the observed 
relic density. Hence, insisting on a DM particle that is motivated independently by particle 
physics naturalness considerations we have delineated a viable scenario that can account for 
the PAMELA observations without requiring a fine-tuning of the parameters, and addressing 
at the same time some of the shortcomings of the MSSM. 
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